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We derive semiclassical expressions for spectra, weighted by matrix elements of a Gaussian ob- 
servable, relevant to a range of molecular and mesoscopic systems. We apply the formalism to the 
particular example of the resonant tunneling diode (RTD) in tilted fields. The RTD is a new exper- 
imental realization of a mesoscopic system exhibiting a transition to chaos. It has generated much 
interest and several different semiclassical theories for the RTD have been proposed recently. 

Our formalism clarifies the relationship between the different approaches and to previous work on 
semiclassical theories of matrix elements. We introduce three possible levels of approximation in the 
application of the stationary phase approximation, depending on typical length scales of oscillations 
of the semiclassical Green's function, relative to the degree of localization of the observable. Different 
types of trajectories (periodic, normal, closed and saddle orbits) are shown to arise from such 
considerations. We propose a new type of trajectory ("minimal orbits") and show they provide the 
best real approximation to the complex saddle points of the stationary phase approximation. 

We test the semiclassical formulae on quantum calculations and experimental data. We success- 
fully treat phenomena beyond standard periodic orbit (PO) theory: "ghost regions" where no real 
PO can be found and regions with contributions from non-isolated POs. We show that the new 
types of trajectories (saddle and minimal orbits) provide accurate results. We discuss a divergence 
of the contribution of saddle orbits, which suggests the existence of bifurcation-type phenomena 
affecting the complex and non-periodic saddle orbits. 



The density of states (DoS): d{E) = 5{E — Ei) of a quantum system -in other words, the density of eigenstates 
Ei at a given energy E- plays a key role in the field of "quantum chaos". Gutzwiller found a semiclassical 
approximation for the oscillatory part of the DoS, in terms of the periodic orbits of the (chaotic) classical system. 
However the DoS, generally, is not probed directly in experiments, as they measure an observable I{E). Often the 
latter can be related to the DoS by a sort of Fermi golden rule: 
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I. INTRODUCTION 




(1) 



In other words, the measured quantity I{E) is the DoS weighted by the expectation value of an observable A over 
the eigenstate with eigenenergy Ei. 
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Experiments probing such a weighted DoS include , among others, spectroscopic studied of atoms in static fields 
1^ , molecular (Franck-Condon) transitions |^ and electronic transport in microcavities . While the Periodic Orbits 
(POs) have offered a powerful tool for understanding these experiments, a quantitative analysis requires one to go 
beyond the Gutzwiller formula. Semiclassical expressions for the "matrix elements" may also be investigated using 
the stationary phase approximation (SPA). 

The first semiclassical approximation for matrix elements involved POs, and in essence reduced to the Gutzwiller 
trace formula (GTF) of the density of states, weighted by the Wigner transform Aw of the observable A, evaluated 
along each PO. This result was found by assuming that Aw is smooth in phase space, and by neglecting it in the SPA. 
On the other hand, the photoabsorption rate of the hydrogen atom was expressed in terms of closed orbits (COs) 
passing through the nucleus because the observable A is very localized. Another semiclassical formula, for the 
conductance fluctuations in microstructures, involved "angle orbits" defined by an angle related to the width of the 
leads [|. 

Hence the the type of contributing trajectories depends strongly on the relative smoothness of A and the semiclas- 
sical Green's function used to represent the DoS. Further, the type of semiclassical expressions also depended on the 
level of approximation used in the SPA integrations. For the semiclassical matrix elements [Q, the variations of Aw in 
both the SPA condition and integrations were neglected. A more refined version developed for molecular transitions 
psf also neglected them in the SPA condition (which yielded periodic orbits), but included them in the integrations. 
Finally, "angle orbits" were obtained by considering both the observable and the semiclassical Green's function when 
applying the SPA. 

A broad range of situations involve Gaussian matrix elements, given by an observable which contains a projector on 



a Gaussian: molecular excitation from a low vibrational state |13 , in the conductance of microcavities with parabolic 
leads |l4| as well as the tunnelling diode experiments which form the subject of the present study. The Gaussian 
matrix elements, conveniently, yield fairly simple integrals, and also have a very clear localization length scale. 

The resonant tunneling diode in tilted fields (RTD) — which is a mesoscopic realization of quantum wells with 
tunneling barriers — was introduced recently as a new probe of quantum chaos . The oscillations of the measured 
current (as a function of the applied voltage) were linked to unstable |^ and stable periodic oscillations, following 
a heuristic application of the Gutzwiller trace formula and taking into account the accessibility of the periodic 
orbits (POs) to the tunneling electrons. 

The RTD experiments generated considerable interest and prompted the development of a series of semiclassical 
theories. One study proposed two separate formulae (obtained by two different level of approximation in the integra- 
tions) using normal orbits starting and finishing perpendicular to the barrier Another study proposed a formula 
using periodic orbits similar to the approach of Zobay and Alber However, it was shown that the normal 
and periodic orbits could not give satisfactory results in all experimental regimes, and that a new type of complex 
and non-periodic trajectories (the saddle orbits) Jl^ provided an accurate semiclassical description. Subsequently, 
Closed orbits, as well as orbits having a minimum transverse momentum were proposed | |TT| in order to achieve similar 
results using only real trajectories. Here we will review and also extend our study of RTD problem, in order to clarify 
previous work and to prescribe the best semiclassical description for this type of problem. 
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In section || we derive a general semiclassical expression for Gaussian matrix elements, expressed in terms of an 
arbitrary type of orbits [eq. (|2^)]. We outline the RTD experiments and show that they are a realization of Gaussian 
matrix elements. 



In section III we explain different reasonable assumptions one can make regarding the localization of the Gaus- 
sian observable in position or phase space. We show that different assumptions yield various types of contributing 
trajectories (POs, COs, normal orbits, etc.), and correspondingly different expressions for the current formula. We 
approach the problem from three levels of approximation: (i) no approximation, yielding saddle orbits (SOs); (m) the 
intermediate formulae, where we neglect one term in the SPA condition; (Hi) the "hard limit" level, where we neglect 
one term in both the SPA condition and the SPA integrations (for POs, this corresponds to the result of Eckhardt et 
al. [0]).Our general formula enable us to reproduce easily any of the five semiclassical theories of the current in the 
RTD that have been proposed in the literature ||9|-[l^. We also propose a new type of trajectories, we term minimal 
orfeitejSelected by requiring that the gradient of the argument of the exponential is minimal (instead of zero as in the 
standard SPA). 

In section |^ we test the different formulae against quantum mechanical calculations and experimental results 
(obtained from Miiller et al. |^ and analyzed in ||2^). We focus on the most interesting regimes, beyond the scope 
of standard PO theory: regimes where there is no real contributing PO ("ghost" regions), or where non- isolated POs 
give non-separable contributions; in these cases the saddle orbits succeed while the PO formula fails. Here we find that 
the closed orbit formula p| represents an improvement over the PO formulae, but requires a complicated strategy 
where one switches between different types of real trajectories in different regimes. On the other hand, the minimal 
orbits achieve the goal of reproducing the accuracy of the complex SOs — across the whole transition from regularity 
to chaos — by using only real dynamics. We also investigate bifurcation-type phenomena for saddle orbits, where their 
contribution peaks sharply, somewhat reminiscent of the effect of bifurcations on the GTF. More detail on this work 
can be found in p^ . 

II. THEORY: DERIVATION OF GENERAL SEMICLASSICAL FORMULA AND THE RTD AS AN 

EXAMPLE OF GAUSSIAN MATRIX ELEMENTS 

A. Semiclassical Gaussian matrix elements 

We wish to find a semiclassical approximation to the quantity 

I{E) = "^{'^|;,\A\tlJ^)5{E - E,) = -^Im Trace [AG] = -^Im j dq j dq'A{q,q')G{q',q) , (2) 

i 

where \ipi) is an eigenstate of the system with eigenenergy Ei. We also introduced the position matrix elements of 
the observable, A{q,q') — {q\A\q'), and of the energy Green's function G{q',q) = {q'\G\q) = lim^^o^ (q'|(i? + 277 — 
H)^^\q). We consider here a closed system with two degrees of freedom q — {x, z), and described by a Hamiltonian 
H. We consider here the special case of Gaussian matrix elements, i.e., projectors on a Gaussian: 

A=m{<f\ , ^x,z)^q^{z)xix) , q^iz)^f^y e-M^^ . (3) 
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can also be a Gaussian (in Franck-Condon transitions) or a cut S{x) (in microcavities or the RTD). To derive 
the semiclassical approximation, one proceeds in two steps. First, one uses the semiclassical expression of the Green's 
function involving all classically allowed trajectories going from q to q' with energy E 



Giq',q) 



E 



q^q' ^ 



iS{q,q-)/n 



(4) 



where S is the action of the trajectories, Px the momentum in a;, and mi2 = Note that we define our monodromy 
matrix M — d{z',p'^)/{dz,dpz) with respect to z and not to the coordinate perpendicular to the trajectory as usual. 
We did not include the phase (— i^7r/2) arising from the number fi of conjugate (or focal) points along the trajectory. 

The second step requires the use of the stationary phase approximation (SPA) in eq. (|^). This states that 
an exponential integral can be approximated by a quadratic expansion of the argument of the exponential around 
the point where the argument is stationary. Depending on the relative smoothness of A and G, different further 
approximations can be made. 



B. Description of the RTD 

A resonant tunneling diode (RTD) can be constructed by adding different layers of semiconductors and applying a 
voltage V between the emitter (where the electrons accumulate before entering the well) and the collector. In effect, 
this will create a wide quantum well between two tunneling barriers (see Fig. |l|) . One also applies a magnetic field B 
at tilt angle 9 in the X ~ Z plane, which creates instability in the classical dynamics in the well. 



conduct ± on 



Jo-^nd ^rof±l& 




FIG. 1. Schematic diagram of the RTD. We show the experimental setup (not to scale), with the conduction band profile 
(effective voltage). Below is the 3-Z) coordinates axis; the magnetic field B is at tilt angle 6 with the electric field —F in the 
X — Z plane. We also show a representation of the distribution of the electrons in the emitter setback; a Gaussian distribution 
(j) m Z due to the magnetic field BcosO, and an Airy function x in X due to the triangular well. The width of the well is 
L = 120 nm = 2267 a.u.. 



4 



The Hamiltonian describing the motion of the electrons in the well can be reduced to two dimensions [p^p[ and 
reads 

1 

H{p,q)^—{j>l+pl)-Fx + —{xsme-zcosef , (5) 
2m 2m 

where we used atomic units (e = mg = 7i = 1), F = V/L and the effective mass of the electron is m = 0.067. The 
length of the well is L = 120 nm = 2267 a. u.. We consider the barriers to be of infinite height; the classical electrons 
will undergo specular bounces [px —Px) at the barriers {x — 0, L), while the quantum wave function ipi of the 
(isolated) well has vanishing boundary conditions [ijji{x — 0, z) = ijji{x = L, z) = 0]. 



C. Bardeen expression for the RTD current 



The current I{E) can be calculated using the assumption of weak tunneling across the emitter barrier, the collector 
playing no important role as all the sites are free for the outgoing electrons (which are accelerated by the voltage 
drop). In that case one can use the Bardeen [|6j formalism, which is a sort of first-order perturbation theory for a 
high barrier 

J^^y\M.fS{E~E^,) , M,^^ [dqU*{q)^{q)-^{q)Mq)]5{x) , q = {x,z) . (6) 



n ^ — ' 2m 7 1 ox ox 

i ^ ' 

In essence, this is an overlap between 'ipi and the "initial state" <&, which is the state of the electron in the emitter 
region prior to tunneling. The overlap is made on a cut taken on the emitter barrier at a; = 0. It was shown p^ ] 
that the use of finite or infinite barriers does not change significantly the important features of the current. Note 
that the Bardeen expression is formally a matrix element like eq. (||), if one writes M.i = with |0) = 

h/{2m)\{pJx + Lpx)®^z\^). 4 = |a:: = 0)(x = 0|. 

In the RTD, one can assume |^ a separable form for $: an Airy function xix) induced by the triangular well, and 
a Landau state (j){z) induced by the effective magnetic field (3 := BcosO. In the experiments under consideration and 
for small 9, one can further assume that only the lowest Landau level [eq. (||)] is occupied 

It is well known that the imaginary part of the Green's function can reproduce the sum of the delta functions of 
the energies, as written in eq. (|^). The current is therefore given by 

J=-^Im j dz j dz' {<^{qW{q')dl,,G{q' ,q) ~ dMnWmdx'^^^^ 

-<i>{q)dx'^*{q')dxG{q\q)^dxm)dx'<^*{q')G{q',q)} . (7) 

Then one can use the semiclassical expression (^ for the Green's function. Because of the Bardeen cut, only trajectories 
going from and to the left barrier q — {x = 0, z) can contribute to the current. Note that the derivatives of -0^ in eq. 
(^ will be transfered to derivatives of the Green's function, yielding factors —px — dS/dx and p',^. = dS/dx'. It has 
been shown [p]Jg[| that the Green's function, as well as its first derivatives, vanish at the hard barrier (this can be 
understood by the Dirichlet conditions for ipi). Hence only the first term in eq. (|^) contributes, and one is left with: 



J=-MEimr3/2x , I=fdzfdz' y f^K^^,^*(^,'ySi-=o,z;x'=o.z')/n (g) 
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D. Derivation of the general semiclassical formula 



First we rewrite eq. (||) using the Gaussian form of the initial state 



dz dz' 



^i^^^')/^' with ip{z,z') ■.= iS{z,z')-^{z^ + z'"^) 



(9) 



-TO12 ... . . 2 

We have also taken the sum over trajectories (a; — Q,z) ^ [x' = 0, z') outside the integrals where it becomes a sum 
over all the different "families" {C\ of trajectories existing in a domain {(z, z') e fi^}. 

The tool used for this kind of integration is the stationary phase approximation (SPA). Briefly, it states that only 
trajectories zq — > z'q having a stationary phase [9(p(zo, z'^^jdz = dip{z(), z'q)/ dz' — 0] contribute to the integrals, all the 
others being either damped by the Gaussian or destroyed by the random cancellations of the oscillations due to the 
action. Then one has to expand the phase quadratically around the contributing type of orbit and integrate, pushing 
the limit of the integration Vl R^. 

We saw in the introduction that a variety of contributing orbits have been proposed in the case of the RTD. 
Therefore, we shall not specify the type of orbits yet, but rather develop a general formula valid for any type, and 



discuss the different choices in section III 



One can relate the second derivatives of the action to the monodromy matrix M , and the quadratic expansion 
of the action reads: 



dS dS 
S{z,z') ~ S{zo,z'g) + Sz — {zo,z'f)) + Sz' — {zo,z'fj) 



2dzdz ^ ^^^ (zo,Zp) + dz -q^{zo,Zq] 



dzdz 



1 



Sq ~ p"Sz + pl,6z' H r^z^mii - 26z5z' + ^'^77192! 

2™-- 



--: S2{5z,5z'\zo,Zq) , 



(10) 

(11) 
(12) 



with 5z — z ~ Z{) and 5z' 



z'q. The "phase" of the initial state is already quadratic. We follow the techniques of 



Bogomolny and Rouben |^ , and complete the square: 



SPA / /3 V- 

(p2h)-= 1^2(7; Zo,Zo 
iSo 
iSo - 



d-y 



-mi2 



^-{z^ + z'^) 



with 



7 = {Sz, Sz') = (z - zo, z' - Zg) 
^ = (-/3zo-zp°,-/3z^ + ip^,) 



. mil ^ 1 \ 



mi2 



mi2 



n 



(13) 

(14) 
(15) 
(16) 

(17) 
(18) 

(19) 



TO12 



.m22 



M2 / 
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£o denotes the different contributing trajectories. With the change of variables 7' = 7 + 7i, the last term in eq. jl^ ) 
gives a pure two-dimensional Gaussian, equal to 27r?i/\/det H. The final result is 



2^ J2 2J^^^^^e[*^°+^(^-^o)]/'' , (20) 



with 

V = -mai - iiSTiM + /J^miz (21) 

B r 

r(z, -^O = - ^ X S 2^ [-'7121 - + [-77121 - «/3m22] + 27/3zz' 

+ 2^Zp^ [7/377722 - ^^77112] + [-7/377711 + /3277712] + 2zp^, - 2z'p^ ^ . (22) 

The above formula describes the oscillatory part of the current. We do not consider here the "smooth" part, 
obtained by considering zero-length trajectories this part varies slowly with the energy, and corresponds to the 
Weyl term in Gutzwiller's theory of the density of states [|I|. 

The loss of coherence due to phonon scattering is not considered in the formalism presented here. It can be modeled 
by adding an exponential factor exp(— T/r) in the sum, where T is the real part of the total time of each trajectory 
and r is the damping (decoherence) time. We shall proceed the other way round, canceling the effects of the damping 
in the experiments. 

The formula is valid only for isolated expansion orbits. Also, we did not write explicitly the phase arising from 
the number /i of conjugate points along the trajectory, as we shall primarily consider individual contributions from 
isolated trajectories. 



III. THEORY: SEMICLASSICAL FORMULAE FOR SPECIFIC TYPES OF TRAJECTORIES 

At this stage one should go back to the SPA applied to eq. (^), and examine which choices of expansion orbits 
have been or can be made. First we mention one remarkable feature of the RTD, which is that for any starting z, 
there exists, generically, a starting momentum for which the trajectory is a time-reversed duplicate of itself and 
therefore closed. We call such trajectories time-symmetric (TS). They are defined by 

TS : (z,p,)^(z',p;) = (z,-p,) (23) 

and satisfy the important relation 77711 — ^22- The existence of TS orbits is a consequence of the fact that for any 
starting z, one can find a starting so that there is either a perpendicular bounce on a wall [pz{x — or L) — 0], or 
a turning point on the energy surface at a point where p = 0. Note that some self-retracing trajectories with ^ 
or z' ^ z are not TS. 

One actually finds that each choice of expansion orbit shown below contains a subset which is time-symmetric (TS), 
and that in almost all cases only the TS subset contribute to the current. This is the reason why we first write (j2^) 
for time-symmetric (TS) orbits. Using (03) and 77111 = '7722, one finds 
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(24) 



where we have defined 



6^-^P^^ . (25) 
mil - 1 

We now consider different possible choices for the expansion points. 

A. Saddle orbits (SOs) 

The first expansion orbits we consider here are given by the strict apphcation of the stationary phase condition on 



SO: <; '"^'^'^ (26) 
We call such trajectories saddle orbits (SOs) [g2[. Inserting eq. ( p6| ) in eq. (|2^), one finds Tsoiza, z'^^) 



.,2 , '2 



2 I ^6 ^0 



Xso= E 27rhJP^e^'''>/'^<t>izo)H4) ■ (27) 

SOa ' 

In the case where SOs are TS (TSSOs) one has rTsso(-2o) — ~P4- We shall show in section |^ that the SOs are the 
most successful types of orbits for a semiclassical description of the quantum current. One difficulty with SOs is the 
fact that they are complex. This is the reason why Bogomolny and Roubcn |9j decided to avoid them (and considered 
real trajectories), although they were well aware of the fact that the stationary phase approximation yields SOs. 

The SOs are non-periodic; this means that one cannot look at repetitions of a "primitive" SO, as it would not 
satisfy the SO condition. Instead one must search for an SO with a higher period. Contrary to complex POs, the 
complex conjugate of an SO is not an SO. 

B. Normal orbits (NOs) 

To obtain real trajectories, one has to make a further approximation and neglect one term in the SPA condition (|2£ 
Bogomolny and Rouben considered that the dynamics in the well are very chaotic; in that case one should expect 
the oscillations due to the action to dominate the Gaussian damping of the initial state. Formally, this corresponds 
to taking the limit /? ^ in (p^), and yields normal orbits (NOs): 

NO : I " (28) 

Essentially, this states that the contributing trajectories are determined solely by the oscillations of the Green's 
function: they must cancel the variations of the action, however small their accessibility to the initial state is. 
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Moreover, TS normal orbits (TSNOs) have Zq = zq- This imphes that TSNOs are actually a special subset of 
periodic orbits, that are time-symmetric (TSPO) and start with Pz=0. In the non-TS case (zq ^ zq), the second 
repetition of the NO is actually a TSPO; the first repetition of a non-TS NO is in a sense a "half PO". One finds 



(zo,Zo) = - — Zq [-m2i - i/3mii] + (zo)^ [-m2i - ?;/3m22] -I- 2i/?zo2;o > , (29) 




NO I 

which is equation (109) of Bogomolny and Rouben(1999) ||]. In the TS case, one has rTSNo(2o) = —jzrgZQ, and the 
current can be written 



^TSNO- 2. 2W e e , 7-^^|^|2 , " /^^o ^ ^ j^p ■ («^U) 



NOs 

(=0-4) 



We shall call this expression the PO/NO formula. Note the shift AS" of the frequency of the oscillations from the 
action S of the PO. Also, the term 7 can reduce the Gaussian damping due to the initial state, for trajectories with 
large zq. In fact, this formula takes into account torus quantization effects a la Miller p9[ | occurring in large islands of 
stability surrounding stable POs. It has been shown [|j that it is analytically equivalent to a model building the 
current as an overlap between the initial Gaussian and the wavefunction in the well approximated by the harmonic 
oscillator state corresponding to Miller tori. 

It is also interesting to consider the case when /3 is very small. The first terms of an expansion of F and D in powers 
of /3 yield 



2^HLNO —1^21 , FhlNO 2 



[z2 + (4)2] ^ Xhlno- V 27r?i,/^e^^/''0(zo)0(4) (31) 

txi \ m2i 



We refer to this kind of expansion as "hard limit" (HL). This was the first formula proposed by Bogomolny and 
Rouben It is justified in the case of extremely chaotic dynamics, where the oscillations of the Green's function 
are supposed to be much stronger than the Gaussian decay of the initial state. In the case of a TS orbit, the precise 
condition for the validity of this theory is [p| 



dz 



mil - 1 



/3mi2 



^»1 . (32) 



The HL result corresponds to the SPA method applied to the Green's function only. As (3 is supposed to be very 
small, the initial state function is neglected in the integrations (|l3|), and is taken out of them; it is evaluated at the 
NO and gives the simple Gaussian factor Phl- The integral is carried out only over the Green's function, bringing in 
a prefactor 77121 and the exponential of the pure action. Note that this is the usual way of proceeding with the SPA 
method, while considering the variation of both functions in the integral ( p^ is not standard. 
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C. Central closed orbits (CCOs) 



Narimanov and Stone proposed a semiclassical approach|^ which effectively amounts to consideration of the 
other extreme case, where the Gaussian damping dominates the action oscillations; this assumption can be justified 
for fairly regular dynamics. This corresponds to taking the limit f] ^ oo in (p6|), and yields central closed orbits 
(CCOs): 

f zo = 

ceo : <^ (33) 
I ^0 = . 

Here the contributing trajectories give maximal accessibility to the initial state, while they do not cancel the variations 
of the action. In this case one finds 

rcco(p°,P°0 = [-if3m22 + f3^m,2] + {pI? [~if3m^i + fS^m^^] +2i/?p°p^,| . (34) 

This formula is equivalent to equation (14) of Narimanov and Stone ||ll|. (They derived a general formula for any 
number of excited Landau levels in the initial state.) TS central closed orbits (TSCCOs) have = — p", and give 

rTscco(^o) = ^T4ib^)'- 

One can also consider the hard limit, i.e. the first order expansion in 0{\/ [3): 



2?HLCCO - /3'mi2 , T^i.cco^^^J{plf + {pl?] ^ ^hlcco - ^ 2nhJ P^e^s/'^4>{p'.)kp'l') ■ 

2/3 cT^, V -mu 

(35) 

Here we introduced the initial state (i.e., the observable) in momentum representation: (j^iPz) = 
(/37r?i)~^/'* exp[— p^/(2/3?i)]. The hard limit is equivalent to neglecting the quadratic term of the action in the in- 
tegral (p^. The integration of the linear term with the initial state is in effect a Fourier transform, and brings in 
the value of the initial state in momentum representation at the CCO. Alternatively, one can express the Green's 
function and the initial state in momentum representation, and argue that the latter is smooth and can be taken out 
of the integral by stationary phase approximation. A similar expression in terms of "closed orbits at the nucleus" and 

— 1/2 

involving a weighting by mj^j "^as found in the semiclassical theory of the photoabsorption spectra of a hydrogen 
atom in external fields The similarity is somewhat limited, as the expression for the photoabsorption spectra is 
much more complicated than mere momentum wave functions (it involves matching the semiclassical Green's function 
to a quantum one in the vicinity of the nucleus). 



^This was a complement to their periodic orbit formula presented in Narimanov et al. and discussed below in subsection 

iiiD. 
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D. Periodic orbits (POs) 



Periodic orbits are a natural choice, as it follows the expansion around POs found in the derivation of the Gutziller 
trace formula for the density of states. A discussion of this choice is more adequately made using a phase space 
formalism, described in appendix Alternatively, a more direct route is to define "average-difference" coordinates 



Z^l{z' + z) \ Pz = \{p'z+Pz 

I\z^ z' - z y Apz = p'^ - Pz 

which one uses to write the z-observable A = \(t>){4>\ in position space as: 



(36) 



-,2 



Aiz, Az) = cj>{z - ^AzW{-z + ^Az) ^ y e--' " , (37) 

while the Wigner transform is defined by 

W{z,pz) = f dAz e'P-'^'/^Aiz, Az) = ^e"^^' " . (38) 

ZTTn J Trn. 

Here we have written two different Gaussian widths (3z and (3/^z for respectively z and Az. Of course in reality 
we have (3z = Paz = /3, but retaining the distinction clarifies the following discussion. The action is S{z, Az) = 
S{z — Az/2, z + A/2), and its quadratic expansion around a point (zo, Azo) reads: 



S2iz,Az) = So + Ap"Jz+fJAz 
1 



6z^{TiM - 2) + Sz6Az{m22 - mn) + -SAz^{TtM + 2) 



(39) 



2toi2 

with Sz — z — zq and 6Az = Az — Azq. 

The idea is to apply the SPA method to the integral I oc / Az / z A{z, Az) exp[iiS2(z, Az)/?i], i.e., with respect to 
the variables (|3^). The SPA condition reads: 

pO = -iz/3A^Azo ^^^^ 
Ap° = -2il3zZo . 

In Eckhardt et al. one assumes the Wigner transform to be smooth as a function of {z^pz). This corresponds 
to the case /Sj — > and /Ja^ oo, which gives for (|40|): 

f Azo = 

POs : <^ " (41) 

that is, periodic orbits (POs). Hence POs arise naturally when one consider a smooth Wigner transform; as 
A(z, Az) is its Fourier transform [see eq. ( |Al[ )], it is smooth in z, but localized in Az. This corresponds to a "local" 
operator, in the sense that A(z, Az) ~ a{z)5{Az). This also enables one to recover the Gutzwiller trace formula, via 
i ^ 1 ^ A(z, Az) 5{Az). 

Things are different for the Gaussian matrix elements, which are written as a projector over a Gaussian state. 
They are the product of two functions depending separately on z and z', and cannot have the property of being 
simultaneously smooth in z and localized in Az: either it is localized in both, or it is smooth in both. One cannot 



11 



change (3z and /?az independently. This fact was noted by Zobay and Alber ||l3[ in their work on Franck-Condon 
molecular transitions, which involved very similar equations. 

Nevertheless, it is still fruitful to consider periodic orbits for the RTD. Putting z' = z,p'^^ = pz in (p^), one finds 

rpo(z",p°) = -^{4 [-2"i2i - ^f3{TTM - 2)] + [-^f3(TrM - 2) + 2f3'm,2] 




(42) 

This formula is equivalent to equation (19) of Narimanov and Stone p^ . An important subset of POs are TS, and 
havep^ = (TSPOs). As mentioned above, TSPOs are identical to TSNOs, and therefore give the same contribution 
®- 

For the hard limit, an expansion in (3z and 1/(3az gives 



I^HLPO-/3(TrM-2)/2z , Fhlpo - - -^(pD' Xhlpo = E tTm -1 "'^^'^^""'^""^ 



(43) 

This corresponds to the formula for semiclassical matrix elements proposed in Eckhardt et al. j7j , that was derived for 
an observable which is smooth in phase space. This formula is basically the Gutzwiller trace formula (GTF) weighted 
by the Wigner transform calculated for each PO.^ To get the hard limit directly from the integration (p^, one neglects 
the quadratic variation of 5 due to Az , and neglects the variations of e^^^^^'^ over the integral (i.e., one uses its value 
at the PO). The integration of e"'^'^^ ^jth the linear term in S due to Az is a Fourier transform, which gives 
the Wigner part e^'^P-^ The integration over the variations of the action due to z gives the TrAf — 2 prefactor, 



as in the GTF. Alternatively, one can work in phase space and apply the SPA to (A3), neglecting the variations of 
the Wigner function W in the integral. Note finally that the hard limit formula for POs expresses in some sense the 
heuristic approach that was first used to interpret semiclassically the current in the RTD, where one considered the 
effects of the stability of POs on the density of states given by the GTF, while taking into account the accessibility 
of the PO to the tunneling electrons 

E. Minimal orbits (MOs) 

Narimanov and Stone [O proposed the CCOs in order to extend the PO formula to regions where one has no 



real POs (also called "ghost regions", see section IV), while avoiding the problem of complex dynamics raised by 
SOs. They also proposed an extension of the CCO formula in terms of time-symmetric orbits which had a minimal 
momentum transfer Apz — —'2pz, i.e., d{Apz)/dz = 0. The argument was that the Wigner transform in the PO 



^The result in Eckhardt et al. contains the average of the Wigner function taken along the path of the PO. In our case, the 
Bardeen cut a,t x — x' = means that we need to evaluate the Wigner function only at the starting and final points of the PO. 
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formula has a Gaussian damping that kills the contribution of trajectories with pz that are not small. This is the 
only proposed semiclassical formula for the RTD that we have not tested by comparison with quantum calculations. 
Instead, we will propose and test another formula which is based on a similar idea. 

The SPA method applied on (^) prescribes finding an expansion point which makes the function ip2{z, z') stationary. 
This can be achieved if one can find a point (zo,Zo) such that the linear term in ([l^) vanishes for all 7 = 
{z — zo,z' — Zq), i.e., ^(zo,^o) ~ (~/5^o ~ iplj—fizQ + ip%) — (0,0). As already mentioned, this requires complex 
trajectories (the SOs). The idea here is to find the real trajectory which minimizes $/^{zo,z'q) and therefore should 
gives the minimal linear term in (|l^). This is in some sense the best real approximation of the complex saddle point. 
Defining 

Ciz,z') = \eaz,z')\' =f3'[z' + z'']+pl+pl, , (44) 

one will look for 

0=i^=^2^_™±Pl (45) 
2 oz TO12 

0='-^ = P^z' + El±P>^ . (46) 
2 oz' mi2 

This prescription defines minimal orbits (MOs): 



. Z' = -7712220 + ^P° 

MO : { ° ^ " (47) 



pO, = f3'^mi2Zo - miip° 

The contribution of MOs to the current will be given by using ( |2^ ) with the {zo,p°; Zg,p^,} of the MO. Again, one 
finds that the most important MOs are TS. TSMOs have p° = zqP^ J^^^^ = zoP\S\. Their contribution wiU be 



mil 

calculated with (B4| 



F. Summary of the formulae 

For the sake of completeness, we mention here the last possibility for neglecting one element in (^). One considers 
the case of a Wigner transform that is very localized in z and pz- This corresponds to /3g ^ 00 and /3az 0, which 



gives for (40): 



. ^0 = 

AO : < (48) 
pPz = . 

We call such trajectories average orbits (AOs), but shall not write nor study their contribution. It is interesting 
to note that for TSAOs, one has zq = z'q = 0, which means that TSAOs are identical to TSCCOs. 

We show in Fig. |2| a schematic representation of the different orbits and their related formulae. We classify 
them according to the level of approximation of the SPA method used in the Gaussian integrations, that is: (i) no 
approximation, which gives the saddle orbits (SOs) and the related formula eq. ([2^); (m) approximation in the SPA 
condition (but none in the integration), which gives the normal [NO, eq. (^o|)], central closed [CCO, eq. (|3^)], periodic 
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[PO, eq. (p2[)] or average (AO) orbits [and their related formulae]; (Hi) approximation in both the SPA condition and 
integrations, which give the hard limit formula (HLPO, HLNO, HLCCO). Then we classify them according to the 
underlying hypotheses regarding the predominance of the Green's function G or the observable A in determining the 
contributing trajectories. This is linked to their relative smoothness in position or phase space. Note that the SOs 
correspond in this classification to the angle orbits found for the conductance of microcavities in the sense that 
both types of orbits are derived without any approximation of the SPA condition. 

NONE i SO 



COND { 



COND 

+ 
INT 



MO 



NO 



PO 



HLNO 



^1 



HLPO 



position 



AO 

ceo 

(HLAO) 

HLCCO 




FIG. 2. Schematic representation of the different semiclassical formulae. The vertical axis describes the three different levels 
of approximation: one can neglect one function {A or G) in the saddle point condition [COND] , also in the integrations [COND 
+ INT], or in none [NONE]. The horizontal plane describes the relative localization/oscillations length scales of A and G (e.g., 
the left end means that the oscillations of G in position space are on a much smaller scale than the localization of A, etc.). 



IV. COMPARISON BETWEEN SEMICLASSICAL RESULTS AND QUANTUM CALCULATIONS AND 

EXPERIMENTS;ANALYSIS OF RESULTS. 



A. Scaling the classical dynamics 



In our comparisons between classical/semiclassical/quantal dynamics we exploit an important property of the RTD 
: its Hamiltonian can be scaled with respect to the magnetic field [|2^. Then, the classical dynamics depends only 
on the ratio e — F/B^ instead of the independent values of F and B (the ratio R ^ E /V is roughly constant in the 
experiments) . The experimental regime j6| corresponds to the interval 1000 < e < 100000. The classical dynamics in 
this range evolve from chaotic (low e) to regular dynamics (high e) [|2^. It is preferable to scale the action not with 
respect to B, but rather with respect the action of the 9 — 0°,B~0,R^1 problem: Sq ~ 2Ly/2mV{R + 1/2). In 
this integrable case, the number of oscillations of the wave function approximated by the WKB method [|l| is given 

by 
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which we shah consider as a measure of the "effective in the general case as weh. We define the scaled action 

by S'(e) := S/Sq. This definition is convenient as the three types of experimental oscillations |^ then correspond to 
trajectories with S* ~ 1, 2 or 3. We term these period-one, period-two and period-three trajectories respectively. Also, 
S depends only on e, but is roughly constant as e varies. We called the important period-one trajectories t and the 
most important primitive period-two trajectories s ||20| . 

We can obtain the period of the voltage oscillations generated by a given trajectory. The frequency of the oscillations 
of the semiclassical current is given by the imaginary part of the argument of the exponential in eq. (|2^): S — 
Re S + Im r. We can also define its scaled version E — J^/2TrJ\f. Then one has two consecutive maxima {V, F -I- AV} 
in the current- voltage trace (neglecting the variation of arg V) when 

27r = + AV") - = A[S] = A [27rE(e)7V(V^)j = 27rS— ^ AV = ^ . (50) 

L J 2V WE 

This can be contrasted to the heuristic interpretation based on the DoS that was used before ||^,^ , where the voltage 
oscillations were directly obtained from the energy oscillations given by the GTF: Ay — AE/R — 27111/ RT, where T 
is the period of the contributing PO. This relation is not exactly correct for two reasons: the current oscillations are 
not given by the action S but by E, and the period T arises in the GTF from dS{E, B, V)/dE = T taken at constant 
B and V, while in our case V varies with E through the constant R. 

B. The different types of orbits 

Because of the decoherence induced by phonon scattering, only the shortest trajectories contribute to the current 
in the RTD we analyzed (which has a width L = 120 nm = 2267 a.u.). We compare in Fig. |^ the shape of the different 
types of contributing (all of the i-type). Examples of plots of the other important class of orbits, the s-type orbits 
may be found elsewhere eg ^2 11|. In (a) we present the traversing periodic orbit (PO) to, which makes one bounce 



on each wall, and is responsible for the broad experimental voltage oscillations |20|j5|]. It is perpendicular on the left 
wall (p° = 0), and is therefore time-symmetric (TS) and also a TS normal orbit (NO). There are two period-two POs, 
born in two pitchfork bifurcations around e = 13000. ty is self-retracing but non-perpendicular; hence it is not TS: 
reversing the momentum at the end of the trajectory (on the left wall), one does not find oneself on the portion of the 
trajectory on which the orbit started, tu is TS and therefore also a TS NO; there is also another non-TS NO hidden: 
it is "half" the PO [zq ^ -80 ^ 2^ = 550) and wiU be denoted by tu-NO. 

As e decreases towards the chaotic regime, disappears with an unstable partner in a tangent bifurcation at 
e — 6500. They are replaced by a pair of POs which are complex conjugates of each other; the one giving a physical 
contribution is called "ghost" and is denoted by igh. At e = 3000 (b), a new real PO ti has appeared. We also show 
the saddle orbit (SO) and minimal orbit (MO) that are related to the t-type trajectories; they do not disappear in 
any bifurcations and arc linked to all the three POs (io,igh) and ti as e decreases. In this region the SO and MO are 
between tgh and ti. 
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In (c) we illustrate the non-periodicity of SOs. We show to— SO, related to the primitive PO to, and the very 
different 2to—SO, that is related to the second repetition of the PO. Because SOs (as well as MOs and CCOs) are 
not periodic, one cannot continue the propagation of a primitive orbit, but one has to look for another orbit with the 
adequate action, period and number of bounces so that it corresponds to the repetition of a primitive PO. Note that 
in some cases (e.g., to~SO a.t 9 — 27°, e < 17000), one cannot find an SO corresponding to the second repetition of a 
PO, although one has the SO linked to the primitive PO. 
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FIG. 3. Shape of different types of orbits of the t-type at = 11°. The right wall is at a:: = L = 120 nm = 2267 a.u.. For 
complex trajectories, we show the real part (Re a;,Re z). (a) e — 14000; the periodic orbit (PO) is self-retracing but not 
time-symmetric (TS), as it not perpendicular on the left wall (p^ 7^ 0); to and tu are TS POs and therefore also TS NOs; the 
"half PO" tu—NO [zo — —80 —> z'q — 550) is a non-TS normal orbit (NO), (b) e — 3000; here we have the complex ghost t^h 
which has appeared in the tangent bifurcation of to—PO at e = 6500; we show the related saddle orbit to— SO and minimal 
orbit to— MO; we also have a real PO ti. (c) e — 18000; we show the saddle orbit to— SO corresponding to the primitive to-PO, 
as well as the saddle orbit 2to—SO corresponding to the second repetition of the PO. (d) e = 10000; the central closed orbit 
to —CCO is defined by zo = 0. 



An example of a central closed orbit (CCO) is shown in (d): it looks very different from the related PO. 



C. 9 = 11° theory and experiments: ghost regions and torus quantization 

The method used to compare semiclassical, quantum and experimental results was explained by Saraga and Monteiro 
. For each e, we generate a scaled quantum (QM) and semiclassical current that oscillates with M , in the range 
12 < A/" < 42 , corresponding to the experimental B range read at constant V — 0.5 V. We Fourier transform the 
current with respect to the pair A/", S and get a power spectrum that has peaks at certain values of the action. 
The height of the peaks gives the amplitude of the oscillation, while their position indicates their type: period-one 
oscillations when 5 ~ 1, period-two when 5 ~ 2, etc. For the experiments, we read the amplitude of the oscillations 
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directly from I — V traces provided by G Boebinger , that we analyzed and presented in |gO| . We correct the 
experimental e by 30% to take into account effect of the voltage dependence of the mass. To allow for damping due to 
phonon scattering, we scale the amplitudes by exp[T/T], where T is the (real part of the) total time of the contributing 
classical orbit and r ~ 0.1 ps |^ is a decoherence time. We use here the value t ~ 0.115 ps estimated by comparing 
the maximal values of the quantum and experimental amplitudes. We normalize all amplitudes to the amplitudes at 
= 0°, where for the semiclassics we have F = and \'D{0 — 0°)| = 25, Ve. 

Fig. ^ presents the amplitudes of the different semiclassical formulae at 9 = 11°. Here we study period-one 
oscillations, which correspond to the broad voltage oscillations seen in the experiments |^ and called "i" series [|]. 
They arise from trajectories making one bounce on each wall; ( the to orbit at high e and ti at low e). 

In (a) we see that the quantum calculation based on the Bardeen model reproduces quantitatively the experimental 
behaviour, over a large range of parameters (corresponding to 3 T < B < 12 T). We did not read experimental 
amplitudes for e > 25000, because of the presence of period-two oscillations. 




1000 10000 1000 10000 

e [a.u.] 

FIG. 4. Amplitude of the semiclassical formulae for period-one oscillations at 6 = 11°, compared with quantum mechanical 
calculations [QM, dotted line with squares] as the dynamical parameter e varies, (a) Comparison between experimental results 
[EXP] and quantum mechanical results [QM]. (b)Quantal results and semiclassical results for the PO/NO formula, to which 
three POs contribute: to, the complex ghost tgh and ti. We also show the contribution of tu—NO, which is not a PO. (c) 
Quanta! results and the semiclassical CCO formula, (d) Quantal results and the semiclassical SO and MO formulae. The figure 
shows that while both the SO and the MO formulae give good agreement over the whole range, the PO/NO/CCO formulae 
give agreement only over a partial range. 

The periodic orbit theory is compared to quantum calculations in (b). For to 7 ^gh and ti we use the PO/NO formula 
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(^), which is the common formula given by POs and NOs. The semiclassical formula is accurate in both the chaotic 
(e < 3000) and regular (e > 10000) regions. In the latter, the semiclassical contribution can be understood by Miller 
torus quantization [T9| . The large stable island of tg supports quantum states that are approximately harmonic 
oscillators (HO) functions in the plane perpendicular to the orbit; this will be discussed in more detail below. We 
also see that the contribution of tu—NO to the NO formula seems unrelated to the quantum behavior. Note the spike 
at e = 6500; this corresponds to the tangent bifurcation where and disappear. It is not a divergence, as the 
complex T) does not vanish. The spike is due to the rapid variation of F near the bifurcation. 

The most interesting region is 3000 < e < 6500, where there is no real contributing PO (the "ghost" region). We 
included the contribution of the complex ghost PO igh, but we see that its contribution is too small [Q. A detailed 
view is shown in Fig. 1(a). There we see that the ghost contribution is too small by roughly a factor three compared 
to the QM results. The saddle orbit (SO), on the other hand, describes accurately the QM amplitudes all across the 
tangent bifurcation, the ghost region and the region where ti~PO takes over [see also Fig. ^ (d)]. Finally, we see in 
Fig. H (a) that the unstable partners of the tangent bifurcations and do not contribute to the current. This is a 
general feature that we observed at any angle 6: only a very small subset of trajectories give a significant semiclassical 
contribution. A study of the amplitudes of period-two oscillation (not shown here) shows that other POs like and 
tu are not related to the QM results although their contribution to the PO/NO formula is significant. 




f»-Po -y — t;-po 



(b) 
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FIG. 5. (a) Details of the period-one amplitudes at S = 11° in the low e region. This is the region where no real PO exist, 
and where the contribution from the complex ghost PO fgh is too small, while the SO contribution is accurate and joins up the 
contribution from the POs to and ti. We also show the contribution of the unstable POs and t^ , which are not seen in the 
QM behaviour, (b) Comparison of experimental voltage periods (line with crosses) with the semiclassical period generated by 
the SO to (solid line). 



The contributions of central closed orbits (CCOs) arc shown in Fig. ^ (c). The main objective of the CCO theory 
as presented by Narimanov and Stone |l^ was to complement their PO theory in the absence of real PO (the "ghost" 
region). It partially succeeds, as its amplitude for 3000 < e < 5000 corresponds to the quantum one. However, it is 
clearly inaccurate for e > 5000 and gets worse in the regular region, where one could have expected the assumptions 
underlying the theory to be respected (in this regular regime, the oscillations of the Green's function should be smooth 
compared to the localization of the initial state). Similarly, the CCO theory is not very accurate in the chaotic region 
(low e). 
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We show in (d) the result of the saddle (SO) and minimal orbit (MO) formulae. There is only one SO and one MO 
corresponding to the three POs tg, igh and ti. Both theories are very accurate and reproduce the quantum amplitudes 
across the whole transition from regular to chaotic dynamics. Actually, the MO contribution is even more precise 
than the SO at very low e. 

Finally, we study in Fig. |^ (b) the frequencies of the oscillations via their voltage period A^. We show the 
semiclassical period calculated with eq. ( |50| ) from the saddle orbit to. We do not show quantum periods, which are 
accurately described by the semiclassics. The theoretical periods underestimate the experimental values by some 10%. 
This however is a confirmation of the fact that to orbits arc indeed linked to the broad experimental oscillations. 

Torus quantization is illustrated in Fig. ^, where we present Wigner and wave functions of quantum states con- 
tributing to the current. The wave functions are approximately separable into a HO state and a WKB-type wave 
function along the trajectory. A/i gives roughly the number of longitudinal oscillations, while the number of perpen- 
dicular oscillations corresponds to the pseudo quantum number k of the HO state. One can use this assumption to 
build a current as the overlap between the initial state and the HO state This is valid for stable POs, and it has 
been shown to be equivalent to the PO |l^ and NO formula in the case of time-symmetric orbits. The Wigner 
distributions show the ring structure associated with HO states. 
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FIG. 6. Quantum state contributing to the current in the torus quantization regime a.t 9 = 11°, e = 15000, labelled by their 
eigenvalue A/i . For the Wigner distributions on the emitter barrier (top row), the vertical axis is z and the horizontal axis 
is —pz', the range is adapted to the size of the (classically allowed) surface of section. For the wave functions (bottom rows) 
the vertical axis is x £ [0, L = 2267] a.u. and the horizontal axis is —z G [—2000,2000] a.u.. For A/",; — 20.894 we also show in 
solid lines the classical structures: to for the wave function, and the main features of the Poincare surface of section (points 
representing trajectories hitting the left wall) for the Wigner distribution. In the first two rows, the torus numbers are (left to 
right): fc = 1,3,2 and 0. 
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The hard Umit formulae for POs and NOs (not shown here; see |Q for a test of the HLNO formula) greatly 
underestimate the contribution of off-center POs, because they cannot take into account torus quantization, which 
gives some accessibility to the PO via outer tori which can extend to the z = region at the center of the initial 
state. Finally, torus quantization effects can explain "jumps" in the experimental current, when the dominant torus 
changes with the magnetic field . 

D. Comparison at ^ = 27°: non-isolated orbits 

We present in Fig. |^ amplitudes of period-two oscillations at 6* = 27° (they were called "peak- doubling" regions 
by Miiller et al. Q as there is a secondary voltage oscillation compared to the broad period-one oscillation). The 
quantum model describes qualitatively the rather broad experimental period-two signal (a). However, it overestimates 
it by some 20%. This would be consistent with uncertainties in the estimate of the decoherence time. r. 



(a) (c) 




e [a.u.] 

FIG. 7. Amplitude of the semiclassical formulae for period-two oscillations a.t 6 = 27° , compared with quantum mechanical 
calculations [QM, dotted line with squares] as the dynamical parameter e varies, (a) Experimental results [EXP]; they have 
been multiplied by exp(r/r), where T is the (real part of the) total time of the contributing classical orbit and r ~ 0.115 ps 
is a decoherence time associated with phonon scattering, (b) PO/NO formula, to which three POs contribute: s', 2to and s\\ 
note the gap 3500 < e < 7700 ("ghost" region) between the take over of s" and the tangent bifurcation where s' appears (s' 
disappears at e = 18000 in a "cusp bifurcation" . (c) CCO formula; note the extension of the semiclassical amplitude down to 
e = 4400 in the ghost region, and the inaccurate contribution of 2to. (d) SO and MO formula; note how both formulae describe 
accurately the quantum amplitude from regular to chaotic across the ghost region (there is only one SO/MO related to both 
POs s' and s"). 
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Two types of orbit contribute to the period- two current, which have roughly double action and period than t-type 
orbits. First, we have the second repetition of t, that is a 2:2 orbit making two bounces on each barrier. Secondly, 
there are orbits of the s-type (1:2), making one bounce on the left wall and two bounces on the right wall, with a 
turning point (where the particle runs out of kinetic energy) in-between. 

First we test the PO/NO formula in (b). The peak of the contribution of 2to corresponds to two successive period- 
doubling pitchfork bifurcations (e = 12600 and 14000) of the primitive PO to- This peak can be understood via torus 
quantization effects: at the bifurcations, the winding angle of the stable Iq reaches the value tt; hence two period-one 
torus series corresponding to two successive k numbers become exactly tt out-of-phase and create a strong period-two 
signal. 

There is a very large contribution from another PO (s') in the same region; it also describes qualitatively well the 
quantum behavior, as does 2to. Both POs give a current with very similar scaled frequencies: S ~ 1.9054 for s' and 
S ~ 1.9154 for 2to. Hence, one should not consider their independent amplitudes, as is done in Fig. |^ (b), but instead 
consider the coherent superposition of their oscillatory current. This will be done below in Fig. ^; however it seems 
that these two competing contributions cannot be easily separated. 

s' appears at e = 18000 in a "cusp bifurcation" [Q, due to the discontinuity of the hard bounce on the left wall 
(increasing e makes the trajectory hit the left wall instead of having a turning point on the energy surface). Then it 
undergoes a synchronous pitchfork bifurcation at e — 13600, with the non self-retracing PO s*; there is no effect on 
the semiclassical current. Finally it disappears at e = 7700 in a tangent bifurcation, below which there is no real PO 
able to explain the quantum and experimental signal until e = 4000, where another PO of the same type (s") gives a 
significant contribution. Hence one has another "ghost" region between e = 4000 and 7700. The low e quantum peak 
(e ~ 2000) is well described by si (a 1:2 PO which has one more cyclotron rotation than s' and s"). Note that the 
large spikes are not divergences. 

The central closed orbit (CCO) formula is shown in Fig. |^ (c). The low e is well described by the CCO si, with 
no spike as the one found in the PO/NO formula. The s' contribution to the high e peak is significantly extended by 
the CCO formula, and follows accurately the quantum amplitude down to e = 4700, where it has an unphysical spike. 
The CCO theory also separates the contribution of s' and 2to, as the latter appears where the former disappears, at 
e = 21000. However, the contribution of 2to is inaccurate. 

Finally, we show the saddle and minimal orbit formulae in (d). For both theories, s' describes accurately the 
quantum results all the way through the ghost region down to e = 4000, where there is no such spike as the one found 
with the CCO. They also separate 2tQ and s' in an accurate way, as 2to describes precisely the quantum results for 
high e too. As for 9 = 11°, the SOs and MOs provide the best semiclassical description. While the success of the SOs 
is expected (as they are the correct saddle point of the SPA integrations), the efficiency of the MOs is again rather 
surprising. Note also the consequences of the non-periodicity of SOs, MOs and CCOs, for which the 2to orbit only 
appears around e ~ 20000, while the to orbits (not shown here) exist at lower e. 

Both POs 2to and s' give important contributions for the peak of the period- two signal at 9 — 27°. To build the 
coherent superposition of their oscillatory current, one needs to take into account the constant phase (— i/i7r/2) given 
by the number /i of conjugate points [shown in Fig. M (a)]. The amplitude of their collective contribution is given by 



21 



the height of the peak of the Fourier transform of the current, and is shown in Fig. @ (b). 
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FIG. 8. (a) Maslov index and (b) coherent superposition ("2fo + s'", solid line with dots, ) of the 2to and s' contributions. 
We also show the contribution of the non self-retracing s, PO, and the QM results, (c) Experimental (line with crosses) and 
semiclassical (solid line) voltage periods. 

The coherent superposition [2fo + s'] is much larger than the quantum amplitudes, because the individual isolated 
contributions are already as high as the QM results. Note the discontinuous change at e = 13600; it is due to the 
discontinuous change of /i for s' at the pitchfork bifurcation. It is clear that the coherent superposition of 2to and s' 
cannot, whatever their relative phases, describe accurately the quantum results. Note that these POs are not involved 
together in a bifurcation; this is not the usual breakdown of semiclassics near a bifurcation, that one could solve with 
the use of normal forms (cubic expansions of the action). Also, 2to and s' seem to be well separated in position space 
(their starting positions are — for s' and zg ~ 600 for 2io)- 

Looking at quantum states contributing to the current reveals that the two POs are, in some sense, not isolated. 
A Wigner distribution and the related wave function arc shown in Fig. 0. 





N=20.065 

FIG. 9. Left: Wigner distribution (color plot) at ^ = 27°, e = 16000 and classical Poincarc surface of section (dots), in the 
{~Pz,z) plane on the emitter wall. The square on lower-right corner represents h. Right: corresponding wave function in the 
{ — z, x) plane, with the POs to and s'. The scales are the same as in Fig. ^. 



The Wigner distribution has the familiar ring structure of a fc = 2 quantized torus in the stable island surrounding 
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to. The ring is nevertheless distorted in some way and is also localized on the (here stable) PO s'. Similarly, both POs 
are within the region of the localization of the wave function in position space. We can conclude that the quantum 
state cannot "distinguish" the two POs, and that the POs are hence non-isolated: they contribute collectively to the 
quantum state and to the current. The use of either SO or MO orbits, however, circumvents this problem and yields 
good results throughout. 

E. Comparison at 9 — 20°: divergence of the saddle orbit formula 

Fig. |o| presents amplitudes for the period-two signal at 6* = 20°. The situation is similar to 6 = 27°; we do not 
present here the CCO formula. The low e quantum peak can be described rather well by the PO (a), SO (b) or MO 
(c) formulae. 

The quantum model can describe qualitatively the shape of the experimental peak at high e (c), over a large range 
of parameters (10000 < e < 40000). The 15% discrepancy is probably due to a small inaccuracy in the estimate 
of the decoherence time r. Semiclassically, we have the same two competing orbits as at 6* = 27°. Although its 
contribution is important, it seems that s' does not influence the quantum amplitudes. Note the difference between 
the s' contribution given by the PO/NO formula from the one given by the SO and MO formulae. 




10000 20000 30000 40000 

e [a.u.] 

FIG. 10. The different semiclassical theories for P2 amplitudes at 9 = 20°. The broad maximum is related to 2to and s' 
orbits, while the lower maximum is given by Si orbits. We show QM results and experimental readings [EXP], (a) PO/NO 
formula, (b) SO formula. The inset is in a larger scale, (c) MO formula and experimental results. 

The contribution of 2to to the PO/NO formula is good, but the position of the peak is not very accurate. The SO 
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formula yields unexpected results: it has a very large peak (see inset), which we shall investigate below. The MO 
formula gives the correct position for the peak, but the height is not as precise as could be expected. 

We investigate in Fig. |ll| the large spike of the amplitude of the contribution of 21^— SO to the saddle orbit formula. 
We see in (a) that the reason for this is the fact that the determinant T) [eq. (pT|)] of the quadratic expansion used in 
the SPA integration almost vanishes around e = 17500— 17800. Both its real (solid line) and imaginary (dashed line) 
part simultaneously approach zero. This is a remarkable coincidence, as I? is a complex function, that one expects to 
vanish only if one can vary a parameter in the complex plane (i.e., two real parameters). In this case, varying e only 
on the real axis approaches very closely the zero of T), which should be reached for a value of e with a small imaginary 
component. 
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FIG. 11. Classical characteristics of the SO 2to around the spike at e ~ 17000 — 18000, ^ = 20°. The solid line indicates 
the real component and the dashed line indicates the imaginary component, (a) Scaled determinant T) = T)/ B [eq. (^)] of the 
quadratic expansion used in the SPA integration, (b) Scaled action S. (c) Starting position zo. (d) Real shape (Re a;. Re of 
2to~SO at e = 17476; for fig. (d) the dashed line indicates the limit of the region allowed by {real dynamics) energy conservation. 



The classical characteristics of 2to —SO in that region are not smooth. The real part of the action (b) reaches 
a maximum value at e = 17588, while the imaginary part changes abruptly over the range 17000 < e < 18000. 
The imaginary part of the starting position (c) behaves similarly. The real part of the starting position reaches the 
minimum value Re zq — —89 a.u. at e = 17476. 

We show in (d) the real shape of 2t(jS0 ate — 17476. The outer leg (with high Re^) hits the left wall perpendicularly 
at a point which is very close (less than 20 a.u.) to the limit surface of the region accessible by trajectories defined by 
real dynamics. This limit is where some self- retracing trajectories (like s' and si) have a turning point (i.e., have zero 
momentum). It is impossible to find trajectories with the same bounce structure (2:2) for Rezo < —100, because they 
would miss the intermediate bounce and go back to the right wall via the limit surface. This is similar to the "cusp 
bifurcation" of the PO s'; it is also the only observed mechanism that can remove an SO, MO or CCO as e changes. 
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As e increases, Re^o decreases until it reaches — 89a.u.. If it evolved smoothly, one would expect the SO at higher 
e to start with a lower Re zq -which seems impossible as we have seen that such starting condition did not allow the 
correct bounce structure. Hence, one would expect "a cusp disappearance" of 2to^SO. This nevertheless does not 
happen: one can still find the SO for higher e, as abruptly increasing Rczq satisfies the SO condition. Hence, it seems 
that we have a "failed" cusp disappearance of 2tQ~S0, precisely in the region where the quadratic expansion becomes 
almost degenerate and where the classical characteristics of the orbit are not smooth. 

This is altogether reminiscent of a bifurcation of periodic orbits, where two POs coalesce as the dynamical parameter 
(e) is varied, and where the quadratic expansion of the action (used e.g. in the GTF) becomes degenerate. However, 
PO bifurcations always involve more than one POs, while in this case we have not observed any neighboring SO that 
could coalesce with 2toSO. 

Note that, strictly, one should not integrate (^3|) on the whole (z, z') plane, but only on the domain ft where one 
does find the proper type of trajectories (2:2). Hence one should cut the integral for z < —100; this would yield Error 
functions. 

This situation raises very interesting questions about saddle orbits. Do they undergo a type of bifurcation? Appar- 
ently no, as a bifurcation should involve several SOs, which we have not seen. What is the origin the quasi- degeneracy? 
It comes from a failed cusp disappearance. Does it has an effect on the semiclassical current? Yes, the current shows 
a strong enhancement. Do quantum results show any sign of it? Apparently no, as the QM results are smooth in 
that region. What techniques can be employed to solve that problem? We tried a cubic expansion of the action, which 
removed most of the enhancement but did not show good agreement with quantum calculations. One could try a 
cut-off in the integral. The first task would be to locate the complex value of e for which V — 0. 



V. CONCLUSION 



The general semiclassical formula (20) that we have developed summarizes in a compact way all the theories that 
have been proposed for the current in the RTD (excluding the changes required by the inclusion of excited Landau 
states and a shift of the Gaussian |l^,|l^). It also shows clearly how the different assumptions on the smoothness of 
the Gaussian affect which type of orbits contribute to the current. We found the types giving the best semiclassical 
description: the complex saddle orbits (SOs) and the real minimal orbits (MOs). It appears that the more standard 
periodic or closed orbits do not succeed in all the situations, in particular in the ghost regions and where one has 
non-isolated contributions. 

The (near) divergence of the saddle orbit formula at = 0° raises interesting questions about the SOs, namely the 
existence of a bifurcation-type phenomenon for trajectories defined by a pair of condition (|2^ ) as restrictive as the one 
defining POs (i.e., giving a discrete number of orbits). The complex dynamics we implemented here are very strongly 
restricted by the definition of the hard bounces on the barriers. One could obviously avoid this point by modeling the 
barriers by soft exponential walls. A preliminary study shows that the "cusp" bifurcations are replaced by standard 
tangent bifurcations; however the search for complex POs and SOs with soft barriers appears to be more difficult and 
again raises interesting questions about the high-dimensional search for complex POs and SOs. 
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The techniques used in this work could easily be applied to other systems involving Gaussian matrix elements, in 
particular Franck-Condon transitions and the conductance of microcavities with parabolic leads (where one further 



assumes that the lowest sub-band is occupied |14|). The importance of saddle or minimal orbits would then depend 
on the localization scale of the Gaussian. 

The authors would like to thank E Bogomolny and D Rouben for invaluable help with the semiclassics and the 
complex dynamics, E Narimanov for communication of results prior to publication, and G Boebinger for unpublished 
experimental data. T S M acknowledges financial support from the EPSRC. D S S was supported by a TMR 
scholarship from the Swiss National Science Foundation. 

APPENDIX A: PHASE-SPACE SEMICLASSICS 

Here we describe briefly a semiclassical formalism in phase space, similar to the one used in and jl^. It is of 
course analytically equivalent to the position space formalism described in section ||. We write A{z, Az) ( ^7|) as the 
inverse Fourier transform 

Aiz, Az) = J dpW{z, p)e-'^^''^ (Al) 
of the Wigner function (pq). We wrote p = pz for convenience of notations. The integral (nq) reads then 



^ ^Y. j '^'^ j '^P j dAzy^^W^(z,p)e^'5.(.-,A.)A-,pA.M . (A2) 
First one can integrate over 5Az by Gaussian quadratures without any stationary approximation, which gives 



I = Y, dz dpJ ^^^W{-z,p)e^^°+a{y,P)yn (A3) 



with 



/- -X -A , ;;-A , ._2TrAf-2 / g _ ^22-^11^ TO12 

p) = -pAzQ + dzAp" + dz — 2\p" ~p + 5z — 

2mi2 V 2mi2 / 



D = 



TrM + 2 

[{TtM + 2)(-pAzo + Ap^Sz) - 2TOi2(5p2 + 26pSz{m22 - mn) + 2Sz^m2i] /{TrM + 2) 
TtM + 2 
8toi2 



where we have defined 6p — p — po. Integrating (A3) with the Gaussian Wigner transform (|3q) yields the same result 



as the integration in position configuration (po[). The phase space formalism might be useful for stationary phase 



consideration, e.g. if one supposes that W is smooth. Neglecting W in the SPA condition applied to eq. (A3), one 
find periodic orbits, and integrating one find the result (|4^). 

APPENDIX B: COMPLEX DYNAMICS 

This appendix presents a discussion of general aspects of complex dynamics in the RTD, which is based on empirical 
observations. It raises interesting questions, such as the definition of hard bounces in complex dynamics, the freedom 
in the evolution, and singularities. 
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We complexify the position q, the momentum p and the time t. We keep the energy E real, as it is a physical 
parameter given by the "reality" of the experiments {E = RV). As we know explicitly the formulae giving the 
evolution q(t),p{t) between two bounces, we simply consider their analytical continuation is the complex plane. 
Computationally, we declare these variables as complex and evolve them following a given complex path for the time. 

We illustrate the choice of the complex time path in Fig. We start on the left barrier with Re x = 0, Im a; = 
and some complex initial condition zo,p^ ed. We evolve the time along the real axis < t < ti, until Re x(th) — L. 
For complex starting conditions, x(tf)) = L + iqb is complex. Then we search for the complex time Tc such that the 
imaginary part of x becomes also zero [Im a; (Tc) = 0]. This is done with a Newton- Raphson algorithm starting at 
to = tf,: 

x{tn) - L 



t 



n+1 



Px {in 



lim t„ 



(Bl) 



We consider that this situation defines a bounce on the right barrier, so we flip the x momentum: 

bounce : 



Rex{T,) = L 
Im x{Tc) = 



Tie — > —Re 



(B2) 



A bounce on the left barrier is obtained the same way, replacing L by above. Then, we evolve again keeping 
Im t — Im Tc, until one finds another barrier; we find the new complex time for which the bounce is "real", flip the 
momentum and carry on. 




(b) 



Imt 



Ret 

FIG. 12. Scheme of the complex classical dynamics, (a) Real shape of a trajectory in the (ReXjUez) plane, (b) Time path 
in the complex plane (Ret, Imt). The dashed lines correspond to the direct time path that goes to the bounce in a straight 



line. The dotted lines represent the Newton-Raphson algorithm (Bl) 



Note that although Re x decreases just before the bounce (t < Tc), Tiep^ > because Re t decreases as well. It is 
important here to emphasize the fact that the trajectory beyond the barrier is only a representation of the Newton- 



Raphson algorithm (Bl), and does not represent any physical behavior. We have never defined what the potential is 
beyond the barrier; hence the particle does not feel any barrier until it reaches the time Tc. This is entirely different 
from the possibility that a particle has to tunnel through a finite height barrier, when it follows an imaginary time 
path or has complex momentum. 

In real dynamics, the functions q(t),p{t) between the bounces are analytical, with no discontinuity, divergence, 
singularity or apparent restriction. Their continuation in the complex plane is therefore holomorphic and single 
valued: q{t),p{t) is uniquely defined, and does not depend on the specific path one has used to reach the time t. 
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Using Cauchy's theorem, one can deform a complex integral along an open path (keeping the starting and final points 
constant), provided one does not meet a pole of the integrated function. Hence, one can choose a "direct" time 
path which goes straight from t — to t — Tc (see Fig. The particle will arrive on the barrier directly with 

Im x{Tc) = 0, and will undergo the bounce without venturing "beyond the barrier". There is clearly an ambiguity 
related to the freedom of choosing a time path: each different time path yields a different looking trajectory [see Fig. 
|l2|(a)]. However, they all have the same global properties such as the total time T, action S, stability matrix M, 
final position and momentum as a consequence of the "analyticity" of the dynamics. As the current formulae or the 
Gutzwiller trace formula (GTF) use only the global properties of a certain type of trajectories, one can employ any 
convenient time path with reasonable safety. Other topics like the study of scarring (the localization of a quantum 
state on a trajectory) depend strongly on the shape of the trajectory, and are therefore more problematic to address. 
We always plot the "direct" time paths; they give trajectories that avoid the region beyond the barrier and are 
more similar to the trajectories found in real dynamics. For instance, a complex orbit looks self-retracing when the 
corresponding real orbit is self-retracing. 



It is our definition of a boimce (B2), and the specification of the number of left and right bounces (the "bounce 
structure") that determine uniquely the global structure of the time path in the complex plane. It has to go through 
the times T^, i — 1,2, .. corresponding to each desired bounce. One could wonder if some points in the complex plane 
correspond to some sort of singularity in position or time (by singularity we mean either a pole or a point where the 
derivative is not defined.) In this case one should not deform the time path beyond this singular time, as this would 
yield a totally different trajectory. For example, Takada et al. have related singularities to tunneling across a 
finite height potential barrier. So far, we never came across such singular points. Once again, we emphasize the fact 



that we do not have a finite height barrier, but define the infinite barrier through the condition (B2). Takada et al. 
[ p3| deformed the time path along the real axis (— oo < t < +oo) into the complex plane; pushing the deformation 
beyond a singular time corresponded to a tunneling event. Hence, going around a singular time changed the structure 
of the trajectory. Our approach here is different: we start at t ~ and must pass through the bounce times T*. 
In a sense, the latter times are "singular" as the momentum changes there discontinuously; however they are not a 
discontinuity in the potential, one cannot go "around" them and there is no apparent barrier to tunnel into. 

In the complex plane, multi-valued functions are behind the singular points where no derivative can be defined; they 
are the points where the different branches or Riemann sheets join. Semiclassical theories involve specific trajectories 
(x ~ 0, z) ^ (x' = 0, z'), which can be "multi- valued" , in the sense that several starting momenta giving several 
different trajectories can nevertheless connect the same z and z' . In this case, the related singular point is a focal 
point where rni2 =0; they appear in the (Re z, Imz; Rez', Imz') plane of starting and ending positions. We observed 
such singularities in contour plots of the action of TS orbits in the starting z = z' plane. As one goes around a 
complex focal point, one connects the different families of closed trajectories and one has a multi-valued action. They 
are entirely different from the hypothetical singularities in the time plane mentioned above. 

It seems that this implementation of the complex dynamics restrict greatly the freedom that they usually offer. For 
instance, we cannot change the bounce structure of a complex orbit by changing the time path defining the evolution. 
A complex trajectory corresponding to a real trajectory having a turning point behaves also as if it has some kind of 
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turning point (it turns back to the right wall without hitting the left wall); one cannot avoid this fact by trying to 
"force" a bounce on the left wall. 
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